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ON  DIGITAL  PROCESSING  VIA  PDQ  FACTORIZATION 


Introduction 


This  report  discusses  a  method  of  compressing  data  bases  that  des¬ 
cribe  large  digital  images.  The  ultimate  objective  of  such  a  procedure 
is  to  effect  savings  in  computer  memory  required  to  process  images  and 
to  decrease  the  amount  of  data  necessary  to  transmit  images  over  a  com¬ 
munication  channel. 

A  digital  picture  (of  photograph  or  video  image)  can  be 
represented  as  a  matrix  consisting  of  rectangular  blocks  of  size  m  x  n 
where  m  and  n  are  the  number  of  pels  (picture  elements)  in  the 
horizontal  and  vertical  directions,  respectively.  To  process  such  a 
picture  most  effectively  (in  terms  of  processing  time  and  quantities  of 
the  pictures)  by  a  computer,  it  is  desirable  that  the  main  memory  of 
the  computer,  during  the  computation,  has  enough  capacity  for  storage 
of  the  entries  of  the  whole  image  matrix.  Unfortunately,  this  is  not 
the  case  for  most  modern  minicomputers.  Therefore,  different  adaptive 
techniques  have  been  devised  to  cope  with  this  situation.  For 
instance,  C.  B.  Moler  and  G.  V.  Stewart  [2]  described  a  matrix 
factorization,  the  so  called  PDQ  factorization  technique  and  applied  it 
to  image  representation.  The  factorization  is  closely  related  to  SVD 
(singular  value  decomposition  technique  developed  by  H.  C.  Andrews  and 
C.  L.  Patterson  [l]).  The  merit  of  PDQ  factorization  over  SVD  is  that 
it  requires,  during  computation,  much  less  arithmetic  and  the  entire 
maxtrix  need  not  be  stored  in  the  main  memory. 

We  have  applied  the  PDQ  factorization  technique  to  various  digi¬ 
tized  pictures,  and  studied  their  approximations  by  different  ranks 
with  fixed  initial  column  vector  and  observed  the  relationships  between 
initial  column  vectors  P.  and  resolution  of  the  approximated  image  with 
fixed  rank. 

In  the  implementation  of  PDQ  factorization,  we  have  developed  a 
computer  program  called  MATRXVAX  written  in  Fortran.  This  program  is 
described  in  some  detail. 


II.  Analysis  of  PDQ  Factorization 

Specifically,  the  PDQ  factorization  aims  to  express  a  matrix  A 
with  m  rows  and  n  columns  as 

A  «  PDQT,  (E) 
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where  P  is  an  m  x  m  orthogonal  matrix,  Q  is  the  transpose  of  Q,  and  D 
is  a  lower  bidiagonal  matrix,  i.e.  whose  only  nonzero  entries  are 

i  =  dj^  if  is  1 ,2,— rain(m.n)  , 

i  s  i  s  2  ,— oin(ra«.1  ,n)  . 

Assume  m  <  n, 


In  our  computer  program,  each  column  vector  of  P  and  Q  is  a 
normalized  one.  Given  a  (nonzero)  initial  column  (or  column  vector) 

from  A  and,  one  can  construct  matrices  P,  D,  and  Q  in  (E)  by  some 
recursive  formulae.  The  important  feature  of  these  recursive  formulae 
is  that,  during  the  computations,  the  columns  a^,  a^»  —  aR  of  A  are 

being  referenced  sequentially.  We  also  note  that,  in  practice,  the 
inexact  arithmetic  due  to  the  roundoff  error  may  disturb  the 
orthogonality  as  well  as  the  normalization.  Note  that  the  choice  of 
the  initial  column  P^  affects  the  whole  structure  of  P,  Q  and  D.  An 

approximation  Ar  of  A  by  specific  P,  Q  and  D  is  formed  as  follows: 


Ar  =  Pr  Dr  V 


where  P.  Q_  are  the  matrices  formed  from  the  first  r  columns  of  P  and 
r  r 

Q  (in  (E)),  and  Dr  is  the  r  x  r  lower  bidlagonal  matrix  with  o. ,  i  a 
1,2,...,r  on  the  diagonal  and  8^ ,  i  =  2,  ....  r  on  the  subdiagonal. 
Clearly  Ap  is  a  reduced  matrix  (of  A)  of  rank  r. 

Moler  and  Stewart  [2]  showed,  using  theoretical  studies  of  the 
Lar.czos  algorithm,  that  Ar  will  be  a  fairly  accurate  approximation  to  A 

even  for  quite  small  values  of  r.  If  one  compares  the  storage  for  A 

and  Ar ,  their  ratio  will  be 
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R(m,  n,  r)  = 


2nr+2r 


clearly  in  order  to  save  storage  it  is  necessarily  that  2nr  <  n  . 


That  is 


2r  <  n. 

Thus  PDQ  decomposition  will  provide  significant  saving  of  storage 
and/or  transmitting  time  if  r  is  chosen  much  smaller  than  n. 


It  seems  that,  in  general,  for  a  matrix  of  low  rank  (its  corres¬ 
ponding  image  is  of  relatively  low  detail),  the  PDQ  factorization  will 
be  effective  in  processing  and  restoring  the  picture.  In  practice,  it 
is  not  clear  how  to  choose  the  value  of  r  in  order  to  have  an  efficient 
representation  of  a  picture  with  less  computer  storage  and 

processing  time.  However,  from  our  experiments,  so  far,  seems  to  be 

a  suitable  threshold  for  the  reduced  rank  r.  Of  course,  the  properties 
that  lead  to  such  low  numerical  rank  remain  to  be  studied. 


Recall  that  [1]  Andrews  defines  a  condition  number  CR(G)  =  d^d^ 

in  the  discussion  of  the  potentially  efficient  representation  of  the 
image  (picture)  in  terms  of  its  eigenimages,  where  G  is  the  original 
picture  (matrix),  d^  is  the  maximal  eigenvalue,  and  d^  the  minimal 

eigenvalue  (|*0)  in  the  SVD  decomposition  for  G. 

Similarly,  we  would  like  to  define  in  PDQ  factorization  a 
condition  number  for  reduced  matrix  of  rank  k:  (^(A)  *  Maxjd^  g,. , 

i  *  1  »  2,  ...  k}  i 

where  d^^  »  max{o1  ,  ...  a  }/min{(a  ,  a 2,  ...,  a^)} 

and  g ^  »  maxJS^  6^,  ...  S^J/minKS^  &2>  . 

A  threshold  for  r  values  would  be  a  value  of  k  such  that  both  (^(a)  and 
Ck+,  (A)  as  well  as  their  difference  are  large  (of  course,  here  the 
quantity  of  largeness  remains  to  be  decided).  Another  way  that  one  may 
measure  the  closeness  between  A  and  Ap  is  by  defining  the  quantity 
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Clearly  if  r  ■  n,  then  D(A,A  )  *  0.  Therefore  a  threshold  that  one  may 

r 

consider  for  the  approximated  value  r  is  to  require  D(A,Ar)  to  be 

sufficiently  small.  We  hope  to  extend  our  efforts  to  include  the  study 
of  the  threshold  of  r  values  in  our  next  report. 

III.  Description  of  the  Matrix  Factorization  Software. 

Because  of  the  limited  memory  of  the  NOVA  800,  our  principal  com¬ 
puting  tool,  factorization  of  the  image  matrix  is  accomplished  by  the 
program  MATRXVAX  running  on  the  VAX  11/780.  A  formatted  image  disk 

file  'MATAP'  is  created  by  the  program  VAXDATA  on  the  NOVA  800  and 

written  onto  a  magnetic  tape  by  the  program  MAGTAP  using  a  blocking 
factor  equal  to  1 .  Program  MATRXVAX  ( 1 )  reads  the  image  magnetic  tape 
file  into  memory  on  the  VAX,  (2)  factors  the  image  matrix  into  the 

T  T 

matrices  P,D  and  Q  of  specified  rank  r  (3)  reconstruct  the  image  PDQ 

T 

onto  a  magnetic  file  and  (4)  writes  the  matrices  P,  D  and  Q  onto  a 
magnetic  tape  file. 

The  image  and  matrix  magnetic  tape  files  created  by  VAX  are 

written  into  disk  files  ’MATAPIN'  and  'MATP' ,  respectively,  by  the 
Program  MAGTAP  on  the  NOVA  800.  The  program  WRTMATAFIN  writes  the 
reconstructed  image  stored  in  'MATAFIN'  on  the  CCMTAL  CRT  and  the 

program  XCMT  writes  the  CCMTAL  disk  file  onto  magnetic  tape.  Images  of 
rank  r^,  where  r^  £  r  are  written  either  on  the  CCMTAL  CRT  by  the 

program  XMTCT  or  Optronics  C-4300  Colorwriter  by  the  Image  Writing 
Software  System  PWREC. 

IV.  Test  Results  (Pictures) 

V.  PDQ  Factorization  for  Blurred  Images 

Denote  the  original  image  by  A  (without  loss  of  generality,  we  may 
assume  A  is  a  square  matrix).  Suppose  that  the  noise  is  modelled  by 
the  presence  of  two  unknown  matrices  N^  and  N^  such  that  for  any  given 

image  A,  B  =  N^A+N^  will  be  the  blurred  image  of  A.  In  order  to  apply 

PDQ  factorization  to  A,  we  first  have  to  determine  the  two  unknown 
matrices.  To  do  this  we  use  identity  matrices  I  and  21  to  define 

B1  s  Njl  ♦  N2  and  B2  s  2^1  ♦  N2. 

Then  B1  and  B2  are  the  two  known  matrices.  Therefore  N^  *  B2  -  B1 
and  N2  s  2B1  -  Bg.  Once  and  N2  are  obtained  and  stored,  then  for 
each  blurred  image  B,  we  apply  (B-N2)  (provided  that  N^-1  exists). 

We  get  the  original  image  A  and  then  proceed  with  PDQ  factorization. 
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RANK  32  ■  RANK  16 


VI.  Concluding  Remarks 

The  above  results  will  be  used  for  further  exploration  of  the  PDQ 
factorization  to  digital  images.  As  part  of  this  work  efficiency  of 
factorization  needs  to  be  improved. 


So  far,  we  have  seen,  through  some  tests  that  the  PDQ 
factorization  did  provide  significant  saving  of  storage  as  well  as 
processing  time  in  the  representation  of  some  ordinary  pictures.  We 
hope  that  further  research  will  shed  some  light  on  the  underlying 
theory.  Our  next  topics  of  studies  will  be  (i)  properties  of  images 
(matrices)  whose  PDQ  factorizations  produce  good  (effective) 
approximations  (i.e.,  the  effective  numerically  reduced  rank  r  is  much 
less  than  the  rank  n  of  the  original  matrix,  say,  rank  r  less  than 
n/4).  (ii)  Applying  PDQ  factorization  to  the  enhancement  of  the 
images. 


VII.  List  of  Programs 


VAXOATA 

08/20/8!  1*139120 

CURRENT  DIR,  VANG 


C  *009  FILE  nameivaxoata 

c  input  oata  to  vax  prop  tape 

PARAMETER  Na 128 , Ma296 
DIMENSION  Ia(N),P(M).I6(M) 
COMMQN/PPaP/AiB 
Data  Aib/0.#1./ 

PFCX)eA*X*0 

ACCEPT  'MATRIX  OIHENSION  •  »,NA 
ACCEPT  'MATRIX  RANK  •  ' $  NR 
CALL  OPEN  12 » 'MATAF* ,2, IER) 

MRlTE(2,!fe)NA,NR 

10  F0RMAT(IX,2(I3,1X» 

ACCEPT  'TAPE  FILE  NUMBER  ■  '  #NFILE 
I2*N A/8 

ICNTae 

1F(NFIIE.EG.0)GO  TO  100 
CALL  SKFLMT  tlER'NPILE) 

100  1LL«0 

CALL  nEaOMT CIER, IA, ILL) 

IF  (IEK,n£,9)GQ  TO  12* 

110  DU  132  J»l»NA 

IbCJ)a2 

130  CONTINUE 

GO  TO  129 

120  CALL  FMLN(lA,IB,NA/2) 

129  00  20  let .12 

M«(I-l)*0*l 
A2«Hl*7 

w«ITE(2,3W) (IB (L) /L*F 1 »A?) 

30  FORMA! (IX, 8(16, 2X)) 

20  CONTINUE 

ICNTaICNTal 

IF (I CnT.Ew.NA) GO  TO  160 
IF (IEN ,E0,9) 60  TO  110 
Gu  TO  1  UP 
160  CONTINUE 

00  9P  lal ,Na 
P(I)*PF(FlOaT  (IJ ) 
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CONTINUE 
Du  00  1*1.12 
*l*(I*l)*B*l 
112**1*7 

-RITE  l8, 70)  (P(L)  ,L*M.*2) 

FORHAT(IX,0(F*.0,1X)) 

CONTINUE 

CALL  CLOSE (2 . If 0) 

STOP 

END 


HAT* XVAX 

00/26/01  10147146 

CURRENT  DIR.  VAN6 


RODS  FILE  NAPE i  NATKXVAX 

hatrix  factorization 

PARAMETER  N*18,NRR*10 

COHHOn/NATR/  OCN.NNR) ,P(N,NRR) »PT (N) ,QT(N) ,PTI (N) ,QTI (N) , 
I ALP** (n) .DETAIN) .JAR (N.N) . IA (N.N) 

OPEN  (UNIT*21. TYPE* 'OLD'. FILE* 'RATA. DAT') 

OPEN f UN1T*23, TYPE* 'NCR « .FILE* 'MATAF, OAT*) 

OPEN (UNIT*24. TYPE* 'NEW', FILE*'NATP, OAT') 

RE AO (21 . 46R) NA| NR 
pufiH*TC2(ia,ixn 
DU  42m  J*1 .256 
DU  414  1*1.42 

m*ci»f)*6*i 

*2**1*7 

READ (21 . 430) (IA(L.J).L*K1»K2) 

FOR* AT (0CI6.2X) ) 

continue 
CONT InuE 
DO  *44  1*1.92 
*1*II*1)*0*1 
*2*A 1*7 

READ (21.4SP) (P(L.1).L*K1.K2) 

FURNAT (4(F6.P.1X)) 

continue 

0tTA(l)*». 

SUH*0. 

DO  2  1*1. NA 

Su«*9uH*P(I,l)*PCI,l) 

CONTINUE 

SuHaSoRT(SUh) 

DO  4  i*t. nA 
P(I.1)*P(I. D/SUM 
CONTINUE 
00  20  1*1. NA 
SuH*0 , 

DO  10  Nil .NA 

SUM*9uM*FL0AT(IA(F.I))*P(K.l) 

continue 

Q(l.l)*sun 

CONTINUE 

Sun *8. 

uU  30  1*1. NA 
SUN*$UN*0(I, 11*0(1.1) 

CONTINUE 

ALPHA (1)*SURT (SUN) 

DO  40  1*1, NA 

U(I»1 1*0(1. l)/ALPHA(l) 


40  CONTINUE 

Jal 

90  JaJ*l 

irCJ.ST.NtObO  TO  230 
00  70  J2*  1  #  NA 
9uH«0» 

DU  00  KllfNt 

SuNa3UNaPLOAT(IA(J2,IO)*QtK#Jan 
90  CONTINUE 

PT(J2)»SU«»4tPH*f4-n*P(J2,J-lJ 
70  CONTINUE 

K2»J-1 

00  100  J2al,NA 
3uHl»0. 

OU  00  K»1,K2 
9UH«0. 

00  60  JJai'NA 
»UH»8UM*PTCJJ)*P(JJ,K) 

60  CONTINUE 

3U"la9UHl+3UM*P(J2*IU 
90  CONTINUE 

PTHJ23  aPT  ( J2)  "SUH1 
100  CONTINUE 
3Uh«9. 

00  110  !■! rNA 
SuHa$UH*PTI (I)*PTl(IJ 
110  CONTINUE 

0ET  A (4) »3UPT C3UM3 
00  120  1*1 • NA 
P(I#JJ»PTI (I)/BETa(J) 

120  CONTINUE 

00  170  J2»l t NA 
3UHa0, 

00  100  K»1.NA 

SUNa9Un*PuOAT(IA(K»J2))«P(N»J) 
160  CONTINUE 

OT  CJ2) aSUN-BETA ( J) *0  C J2, J-l ) 

170  CONTINUE 

K2»J-l 

00  204  J2»l i NA 

SUN  1  a.), 

00  100  A«1 * A2 
SuHa0( 

00  130  JJal.NA 
SUNaSuHaQT  (JJ)«0(JJ»K) 

100  CONTINUE 

SU" 1 *8UN 1 ♦3UN*0 ( J2 » K) 

190  CONTINUE 

OTI (J2) aQT  CJ2)»SUH1 
200  CONTINUE 

OU  210  t«l »NA 
9UMa9ON4OT2tI}a0T2(2) 

210  CONTINUE 

A|.PNA(J)a306T(9UH) 

00  220  !•! »NA 

QCI/ J)aflTI(I) /ALPHA (J) 

220  CONTINUE 

00  TO  90 
230  CONTINUE 

DO  336  Ml,NA 
00  329  lal.NA 
3UH«0. 

00  310  Jal.NP 

3uHa3UH*P (I #  J) • ALPHA CJ  J  *Q  (K ,  J) 
IP(J.EQ.1)G0  TO  310 
9UH09UH«P (2 , J) *6ET A (J) *Q (ft , JaD 
310  CONTINUE 

I AN (I*K)*IP2X (SUM) 

300  CONTINUE 

00  330  Jtat,290 


* 


00  <130  11*1,32 

«l*lXl»l)*04l 

n2*m+7 

4RITE(23,43l)  <lAR(L,Ji)  ,L*Kl.«2) 

431  FORMAT (1X,8(X9,2X)) 

330  CONTINUE 
I  ACCEPT  340,  JVN 

R£AO (11 #34M) JYM 

340  FORMAT ( 'CONTINUE?  (i*YES*0*NC)  M3) 
IF(JVN,EG.0)GO  TO  300 
X  ACCEPT  390. N0I 

REAu(ll,330)NRl 

330  FORMAT ( 'MATRIX  RANK  •  M3) 

300  CONTINUE 
X2*nR/4 
DO  370  1*1, 12 
Kl*CI-l)*4*l 
»2*w 1*3 

WRITE  (24,300)  (ALPHA(L) ,L*Ki,K2) 
WNITEC24.380) (BETaCLL) ,LL*K1#*2) 

300  FORMAT (1X,4(E12,S,1X)) 

370  CONTINUE 

NH2*NH/4 
00  400  jai,290 
00  400  1*1, NR2 
M«(I-1)*441 
K2**i*3 

WRITE (24,300) (P(J, L) , L*K1 ,K2) 

WHITE (24,300) (Q(J,LL),LL*K1,K2) 

400  CONTINUE 

IF(JYN.feU,M)60  TO  390 

NN*NP ) 

GO  TO  230 
390  CONTINUE 

x  close (unit*2i ,oispc3E* 'qelete  • ) 

X  CLOSE (UnIT *23, DISPOSE* 'SAVE  1 ) 

x  close (unit*24, dispose*  'save  1 ) 

STOP 

EN0 


wrthatafxn 

00/20/01  10140117 

CURRENT  01R,  TANG 


C  ROOS  file  names  urtmatafin 

C  "RITES  image  from  VAX 

C  NOTES  nhEN  FORMAT  written  WITH  (IX,  )  TO  VAX 

C  WORO  MUST  BE  REAO  WITH  NOVA  USING  (IX,  )  FORMAT 

C  "ORu  MUST  SE  mEAO  WITH  VAX  USING  (  )  FORMAT 

C  I.E.  iX.  UELETEO  FROM  FORMAT  STATEMENT 
Parameter  M*313,M*a36 

DIMENSION  I A  (N) , IB (M) 

CALL  0PEN(1, (MATAFIN* ,2,IER) 

CALL  UPCNCT(XER) 

ACCEPT  '  Tu  INAGE  0  •  MOUT 
ACCEPT  'OISPLAY  CODE  •  MCOOE 
accept  'matrix  dimension  -  ',na 

CALL  0SPL(XER,IC00E) 

Ou  20  J*1 ,230 
OU  10  1*1,32 
*l*(2-l)*0*l 
X2*K1"7 


9 


nnn 


•  KAO (1,30)  CIA(L).LaKl,K2) 

X  *»ITE(W,J0)  (IA(LL)#LLaKl/K2) 

30  PORMATClX,0(I0,2X}) 

Ou  «tf  MaK1.K2 
1P{IA(m).lT,0)IA(H)b0 
IFCIA(M),6T (299) IA  (M) »253 
40  CONTINUE 

10  CONTINUE 

CALL  *>ACft(lA'I0,NA) 

CALL  HTINABE CIER, I0UT, J, XB.-NA/2) 
20  CONTINUE 

CALL  CLOSE  Cl, IER) 

STOP 

ENO 


MATPQTP 

00/20/01  10140130 

CURRENT  OIR,  VANS 


ROUS  PILE  NANElNATPOTP 

INARE  CONVERSION  SV  MULTIPLICATION  OP  P,0  ANO  0 
MATRIX  RANK  NR0  RUN  ON  VAX  MUST  0E  MULTIPLE  OP  4 
PARAMETER  N«1024.M«2040|NN*290>NMal2S 
OIMENSION  P(N) ,Q(N) , ALPMA(NN) ,8ETA(NN) ,IA(NN) , 10 CNN) 
COMNON/XNTSl/IP(M) 

COMMON/ XNTS2/ 10  (M) 

euuzvalence  ciPcn.pan.ciotD.Qtin 

ACCEPT  «NEN  RUNT  (laTESf 0aNO)  •  t,JTN 
ACCEPT  'MAX  MATRIX  RANK  NR0  a  *,NR0 
ACCEPT  'MATRIX  RANK  •  '/NR 
ACCEPT  'TO  MTSI'/IPL 
CALL  0PENC3, »MATP',2,IER) 
nlInES*1»24/NR0 
NTRAN0NNM0/4 
*008*1024 
N0LKSB0 
I2*NR0/4 
00  30  1*1/12 
K|* ll*l J *4*1 
K2»Kla3 

REAOC3/20)  (ALPHA (L)  ,LaKl/ft2) 
REAO(J/20)(0ETA(L)/LaKl/X2) 

20  PORnAT(XX/4(E13,0,1X}) 

30  CONTINUE 

IP(JVN, £0,0)00  TO  32 

CALL  CFIL*('PMAT',3,200,IER) 

IP (XER,nC,1)TTPE  'PILE  ERROR* 

CALL  CPILM(IQMAT'/3,2R0/IER) 

IP CIEN.nC, 1)TVPE  'PILE  ERROR  2* 

32  CALL  0PEN(l/'PNAT'/2/IER) 

CALL  OPEN (2, »QMAT»,2,XER) 

IP(JVn.EQ.0)GO  TO  02 

XMAXaNNf>S/4 

00  SP  J*1 /NTRANS 

00  40  1*1 « IMAX 

K1*(I»1)*4*1 

K2*K1«3 

REAOO/20)  (P (L)  / L»K1  *K2) 

REAOC3/20) (fi(LL) »LU«Kl#K2) 

40  CONTINUE 

CALL  "N0Lft(l/ (J-1)*NBLKS,IP,N8LKS»:ER) 

CALL  "RBLK (2/  (J-l)«N0LK3.!Q«NBLKS/l£R) 

X  TYPE  'IP(X) /I*l«32  • i (XP(IXX) « IXX*t«32) 

X  TtPE  'IO(I), 1*1,32'# (I0(IXXX),IXXXal,32) 


10 


31 

92 


209 


140 


150 


160 


CONTINUE 

continue 

CALL  RLSf  C'NT0»,ICf»> 

CALL  ,XE0) 

If <IPL.fcC.0)GO  TO  200 
CALL  S*FLhT(ie»,irU 
CONTINUE 

DO  160  a»1»nTRAN» 

CALL  «O0LK(2, (K-t)*N0LN9,ZQ>N0LKa,XER} 

DU  160  Al»l ,NLXNCS 

KJ»l*l*l)*Nk0 

Xi«L 

OU  IS*  I ■! , Nl  R AnS 

CALL  kOBLNCl, (I»t)«NBLKSf I0,NBLK»,ie»O 
DO  15«  1 1  •  1  # NLlNtS 
I J*  (1 1*1 1  *NRB 
SUH»P. 

DO  1 40  j»i,Nk 

SUN»SUM*P  IX J*J)« ALPHA tJ)*U(KJ*J) 

IHJ.e0.1160  TO  140 

9UM»SU**P l I J*J) *BETA (J)*Q(KJ*J*1) 

CUNT InUE 

Ia(11)«IPXA(SUH) 
IP(lAlU).LT,01IA(Xn>0 
IP(lACXl) . CT ,255) I A (III *299 
CONTXNOt 

CALL  PACK II  a,  XBi 256) 

CALL  *PXTnTCIEH,IBi-126) 

continue 

CALL  rTENOF(IER) 

CALL  CL06E (it IEP) 

CALL  CLG3t(2,IER) 

CALL  ULUSE  C J» IER) 

CALL  «L5El'HT0» ,XER) 
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